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SUMMARY 

Closed-form equations for the lift, drag, and pitching-moment coeffi- 
cients iif two-dimensiona 1 airfoil sections in steady subsonic flow are 
ohta’neu from published theoretical and experimental results. A turbulent 
boundary laver is assumed to exist on the uirfot 1 surfaces. The et foots of 
section angle of attack, Mach number, Reynolds number, and the specific air- 
foil tvpe are considered. The equations are applicable through an e.ngle-of- 
attaek range of -ISO® to +180°; however, above about '20°, the section char- 
acteristic . are assumed to he functions onlv of angle of attack. A computer 
program is presented which evaluates the equations for a range of Mach numbe 
and angles of attack. Calculated results for the NACA 23012 airfoil section 
are compared with experimental data. 


NOMENCLATURE 


c airfoil chord, m 

c, drag coefficient 

d 

Cp friction drag coefficient 

0 ^ mean skin friction coefficient 

c lift coef f ic lent 

c maximum lift coefficient (first lift-curve peak,' 

max 

c m pitching-moment coefficient about quarter-chord 

c m ^ pitching-moment coefficient at a * 0 

c s pressure (form) drag coefficient 

correlation parameter for airfoil drag due to lift 
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I, .lit toil sect ion jut (motif, m 

M M.u'ii tuimbor 

M j-j dr.iv’, divergence Mu'll number 

M M1) mor.ii n\ I ivori’.onv'o M.ioti number 

M , K*,uh number at which t rend of i :n- re. is i r.g lift-curve slope with increas- 

ia>; M reverses; also, high'st Mach number lor which . is 
invariant with M 

Koynolds numb or 

R., Rov s'lds mimbor bast’d on tree-stream velocity 

moan valao of airfoil pressure tool I loiont 
t alt foil soctfon maximum thickness, m 

i airfoil soetion angle of attaok, dog 

a, cutoff a, abovt which the slope ot cj with M Is assumed to bo 

constant, dog 

a a correspond Ing to highest value of H dog 

DI) 

a u corresponding to highest value of deg 

MD 

» reference angle about which dcj/JM vs a curve Is assumed to be a 

mirror Image, deg 

1 « 

'st iLl a<rfoi l section stall angle, approximated as jjT'/j’, *" -\i» 


a^ angle of attack for zero lift, deg 

t ratio of specific heats 


INTRODUCTION 


Piloted simulation of helicopters or other rotorcraft In real time Is 
paced by the speed of solution of the rotor dynamic equations of motion. One 
method of increasing computation sp?ed is the use of a hvbrld (analog plus 
digital) computer, with all totor calculations being accomplished in the 
analog computer. Unfortunately, this may preclude the normal practice of pro- 
viding rotot blade airfoil section characteristics in tabular form, since the 
time required for the very large number of table lookups may negate the 
expected sDeed advantage. Therefore, it was determined to derive closed-form 
equations .or the lift, drag, and pitch!, .g moment of airfoil sections which 

2 


I' 



would exhibit, insofar as possible and practical, the effects of angle of 
attack, Mach number, Reynolds number, and the specific section type. Pub- 
lished airfoil data were used to establish the form of equations which both 
reflect general trends and accommodate specific section type characteristics 
where possible. The coefficients of the equations must be obtained from 
experimental data for the particular airfoil of interest, but the data 
reauirements are not as extensive aa for construction of airfoil data tables 
typically used in helicopter rotor analysis programs. 

The equations presented herein are equally applicable to fixed-wing prob- 
lems, and should be useful in an- case where it is either impractical or 
impossible to use an airfoil se:tion dita table. Conversely, if a particular 
data table is desired but not available, the equations provide a means of 
generating the table (with angle of attack and Mach number as parameters) from 
limited input data. 


ASSUMPTIONS AND LIMITATIONS 


The equations presented herein were derived from experimental airfoil 
section data obtained under conditions of two-dimensional, steady-state flow. 
The airfoil surface condition is assmed to be representat iv^ of in-service 
conditions; that is, smooth, but with a turbulent boundary layer over essen- 
tially the entire airfoil. Considering the combined effects of manufacturing 
defects, service wear, environmental deposits, high Reynolds number and f ree- 
stream turbulence, it is reasonable, if slightly conservative, to assume a 
fully turbulent boundary layer for all full-scale rotors and wings for which 
no special laminar flow apparatus is provided (ref. 1). The assumption of a 
turbulent boundary layer establishes the form of the incompressible drag 
coefficient (appendix B). 

The subject equations are applicable only for flight conditions in which 
the free-stream Mach number is less than one. The Reynolds number, based on 
section chord, is assumed to be greater than 1<10 6 . Certain experimental data 
for the particular airfoil type under consideration must be available. Two 
examples of such data are: (1) moment coefficient at low Mach number and zero 

angle of attack; and (2) at least two values of drag divergence Mach number as 
a function of section angle of attack. All such requirements are presented in 
a subsequent section. 

Equations are provided in this report for section lift, drag, and pitch- 
ing moment coefficient through an angle of attack range of il80 o . However, 
above about +20°, these coefficients are functions only of angle of attack, 
and are independent of Mach number, Reynolds number, anJ the airfoil type. 

This simplification is a consequence of the very small amount of experimental 
data available at high angle of attack. 

It is assumed that the equations presented in this report will be applied 
only to airfoils whose characteristics follow the same general trends as the 
experimental data used to derive the equations (see the appendixes). For 
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example, the drag divergence Mach number should bo .it least approx i mat o 1 
linear witn section angle of nttaek. 


AIKV011. SKCTION KOKl'F AND M.IMI'NT 1 OI'AMON 0 


Kquations for < In* lift, drag, .uul p i t e !i i tip moment cl airloi! sc, t ii‘ 
thionqh •180" angle of .it took, derived in append i xi". \. I', .mil reaper 
arc piesonted bo low. The oipi.it ions .iro pr i no t pal 1 v omplrlo.il in nature, 
in tlio low .untie of attack range reipi i re as input experimental ilata ter 
part lc liar airfoil tvpo. These requirements are discussed in a subs, quo 
section. Assumptions and limitations are given In t he preceding section 
facilitate reference to the appropriate section of the appendixes, equal 
numbers from the appendixes are retained here. 


- 


Lift Coefficient 

Small angle of attack (u i * stt ]p - 


dc . 

l 'v ■ da 1 (1 - V ' 


'stall 1 


(«)-(,') {' * '• Ale [■' " ° 


+ 0 . b ( i. - n 


]) 


where 


A - M’ 


i ! 




If M ' Mj. use 


* '* <£)J ♦ r+b; [*<» - 


]} 


1) + 0. 6 ( ti ‘ - !)• I - (0.4M(M - M 


In lieu of equation (A2), but set dc^/da * 0 OS as a lower limit. 
“0 “ («o^l • 0 i M < M, 


a 0 “ ^‘*0^1 


(-„) 1 “ ('’(l ) r 


_ (M, - M) , M > M, 


1 


max 


C, + C r M 8 CjM*’ + C (( M 3 + C,m‘* + (C ti + c: 7 M C8 )sin(C<, + C^MT 

4 


ns 

t i v i 1 v . 

and 
t lie 
at 

. To 
i on 


(A1 > 


(AJ> 


(AT) 


(A4) 


(A5) 

(A6) 


a 


i 



For negative angle of attack, the form of c c is assumed to he given by 

‘max 

equation (A6). The sign of Co will be negative, and for cambered sec- 

tions the constant term Cj -t i L 1 have a different magnitude also. 

Large angle of attack- For positive angles of attack, 

e - 0.811 

' max 

c, - 0.81 . . „vr— - "> • Valf ' ■’ 2 ” ■ <A1 

vS t a 1 1 


max 

'stal 1 dc . /d i + 
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c t - -1.1 + 1.78(0. 01745(a) - 2 . 156 ) ^ , 90° L , 160° 
Cj = -0.765 , 160° < a - 172.5° 
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For negative angles of attack: 
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i'rag Coefficient 
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For angles of attack less (more negative) than the reference angle n R , 
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in equation (B8). Also the negative angle ct which dc^/dM becomes constant 
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bitching Moment Coefficient 
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INPUT DATA REQUIREMENTS 


The input data that are necessary for evaluation of the lift, drag, and 
pitching moment equations can be obtained from airfoil section test results. 
The specific input variables, and the equations in whicu they are used, are 
listed in table 1. (It should be noted chat S A , L/c , e s /op, and K, which are 
used in the drag equations (B3) and (B5), can be obtained from figs. 7-10, 
appendix B, respectively, ar functions of section thickness t/c.) 

Minimum experimental data (for a particular airfoil) necessarv to evnl— 


uate the Inputs listed 

in 

table 1 are: 


c vs 

a 

for several 

M 

C d vs 

M 

for several 

a 

c m vs 

M 

for several 

a 


The effort requi-ed for data preparation will be reduced if the following 
curves are also available: 


vs M 


max 


M dd vs a 




c vs a, for a low M 
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The c,< vs a curves should extend to at least the stall angle. No test data 
past stall are required because the airfoil section charact or f st ies at . i *’ 
angles ot attack are assumed to be independent of sec ti van tvpe. I he appen- 
dixes provide a guide to t ho use of the above data in tea* prej'iratb'u of i u a . 


EXAMPLE CASE 


The equations presented above were used to calculate lift, drag, and 
pitching moment coefficients for the NACA 23012 airfoil section. Input data 
were obtained from wind tunnel test results for the 23012 given in refer- 
ences 2-4, and are listed in table 2. The calculations were done by a digital 
computer program written to evaluate the equations over a range of Mach number 
and angle of attack. A FORTRAN listing of the program is presented In 
appendix 0. Section coefficients were calculated through an angle of attack 
range of -1P0 C to +180°, and through a Mach number range of 0.0 to 0.9. From 
these results, c , c , tnd c at M = 0.1, tor c, = —20 to +180 , are plotted 
in figure 1. ' d 

Calculated and measured (ref. 2) section aerodvnamio eba rac t e r i s t ic s at 
low angle of attack for three Mach numbers, are compared in figure 2. Overall 
good agreement was realized, although some details of the 23013 behavior were 
not rep rod in cd , since the equations are intended to be sufficiently general to 
represent most airfoil sections. No effort was nude to tailor the inputs to 
improve the correlation. Experimental data for the 21012 .it high ingle of 
attack were not available. Figure 3 presents measured and calculated section 
coefficients as a function of Mac! number. in this case, the experimental 
data are from reference 4. Again, generally good agreement was found, 

CONCEIT) INC REMARKS 


The airfoil section equations presented In this report were developed 
principally by fitting curves to a relatively limited set of experimental 
data. Some ac uracy was sacrificed to obtain general applicability to a range 
of airfoil types. The utility of such equations can be Judged onlv in the 
context of their intended use. If high accuracy is not required, these equa- 
tions provldr the ability to calculate airfoil section aerodynamic coeffi- 
cients at anv angle of attack, over a wide range of subsonic Mach numbers. 

The principal application of the equations will be In fixed or rotarv 
wing computations for which it is not practical to use graphical or tabulated 
section data, or for which a data table must be generated from limited input. 
The equations also provide a method for calculating section lift, drag, and 
pitching moment at verv high angles of attack (t.e., -180° < a ^ 180°1. 
Experimental data of this type are available for very few airfoils. Another 
important application is in the prediction of section characteristics at high 
subsonic Mach numbers in cases where such data are incomplete. 
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AIMM'NI) I X A 


PK' ] VAT I ON OK 1.1 KT (R'KKKKll KNT FOK.UIONS 
S.-IAi 1 ANCl.K OK ATTACK 


In flu* section .ingle of attack i.inge 0 •_ ,. » ! i 22 , ellects n: comp 
sibilitv, Reynolds number, and the specific airfoil section of interest a 
considered, Kor angles of attack below the stall, the section lilt ■ nrve 
assumed to he linear: 


res- 
re 
i s 
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stall 


(A! 1 


I.itt-Curve Slope 


Mach juimber and tji l c k nes_s_ efticts- Kaplan's role y r t . 11 is used to 

represent the effects of compressibility and airteii section thick'.u-s on !i;t 
curve slope: 
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where 



-- = section thickness to chord ratio 
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> * 1.4 for air 
so that 




t/c 

1 + t/c 



1) + 0. 6 (;r 



(A2) 


The utility of this relation was investigated by comparing calculated and mea- 
sured (ref. 2) lift-curve slopes for several airfoils over a range of Mach 
numbers, as shown in table 3. Lift-curve slope at M above 0.3 was obtained 
from the measured value at M = 0.3 and the ratio 


(dc^/da) 

(dc^/da) 


M>0. 3 
M=*0. 3 
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ra 1 i-u l a t ed with eqn.it ion (A2). The di! errnce hctvivn t h-' i . i ! • I an 1 
measur'd 1 i* t- curve slope was 'll' or le-n ,• oept for or o '• < i •: top 

airfoil. 

For the airfoils considered in table i. rnd other a tfoils pr ••■’''.ted in 
reference 2, the trend of increasing ift-’urve slope with 'Is eh nicer 
reverses above approttitrui' el v M - 0.8. Fewever, there is no d i s< ernihle p it - 
tern to the rate of doe 1 i no . References 3 and 6 present "synthesized" lirfoil 
section data (for the NACA 0012 and 0015, respect ively » which do exhibit a 
smooth dee'i itie of litt— carve slope above the levers, i . This snootnicss is due 
to the manner in wlitc' 1 the airtoil characterise ics we'; derived: an iteration 

between successive assumption of airfoil data values , nd comparison oi calcu- 
lated with measured rotor performance. Tin average ate ot decline of lilt- 
curve slope for the two sets ot sect ion data is 

' (Jo , l A i) 

. ^ — = -0.45 per Jo: 

Let M at which the trend of 'the lift-curve slope reverses be M. . Then, 
for M above M. , the lift-curve slope will b> 
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i + t / c 



per deg 


n + o . b ( 



(I S'! (M - M, ) 

(All 


It is necessary to set a lower limit for equation (All. Again, the data avail- 
able (o.g., ref. 2) are not sufficiently regular to provide a trend. There- 
fore, the lift-curve slope for the highest Mach number presented was measured 
for several airfoils in reference 2. The average vale- is 0.05. Thus the 
lower limit of equ*rion (A3) is assumed to he 

dc 4 

2 0.05 per deg 

Reynolds number effects- Airfoil section lift-curve slope Is a weak func- 
tion of Rjc,’ up to 1 > 1 tr or 2'1G' , and is essentially independent of Rjj for 
higher values (refs. 2 and 3). Therefore, for most uses, lift-curve slope can 
be assumed to be constant with Rjg. The value of incompressible lift-curve 
slope used in equation (A3) should be selected with consideration for the 
likely range to be encountered. This value is also a function of the 

specific airfoii section considered. 


Angle of ?ero Lift 

Mach nu mber effects - The variation of airfoil section angle of zero lift, 
«0» small with M and mav be neglected, until M reaches a high subsonic 
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value . Then, for most airfoils, u 0 decreases in magnitude with further 
increases In M. "able A presents measured u 0 values for M above M,, 
v re M, is the highest M in the experimental data for which u,, remains 
c t a constant value. (The Mach number' Mj is also the Mach number at which 
the trend of increasing lift-curve slope with increasing M reverses.) The 
data in tnbie A provide no consistent trend. For six of the sections, there 
is a strong decrease in the magnified, of ii 0 with M, for M above Mj. 
However, for three of the sections, the decline of ii i; is slight, weile for 
two other sections, the decline of u,, is reversed as M continues to 
increase. , Since a general expression for u 0 was desired, it was assumed 
that ttie variation of section angle of zero lir't with M can be adequately 
represented by a straight line above M ( : 

it j = (.ip ) j , 0 i M < Mj (AA) 


a n " <V» “ 




( a n > ? 


Mi - M ? 


~ (M, - M) , 


M > M. 


1 


(A5) 


Mach number M., is some convenient M % Mj for which i () can be determined 
from the experimental section data. Note that for symmetrical airfoil sec- 
tions, * 0 even for high Mach numbers. Thus, for symmetrical sections, 

■\) = (-oil - < rt n>.; = °- 

Reynolds number effects- The angle of zero lift (low Mach number) of an 
airfoil section is determined by the camber. The extensive section data pre- 
sented in reference 3 indicate that is not significantly affected bv 

Reynolds number. 


Maximum Lift Coefficient 

M ach number ef fect s- In the angle of attack range bflow the first lift 

peak, the variation of airfoil section maximum lift coefficient, cr , with 

max 

M may be categorized as: (1) throughout the range of interest, c 

'max 

decreases with M, or 1 2 ) in part of the range, c 5 increases with M 

l mnx 

(but decreases otherwi.se). Experimental data for both types are shown by the 
solid curves in figure A, *nd both can be fitted by polynomials of the form 

c, * C, + C,M + C^M 7 + C,M 3 + C.M 4 + (C, + C 7 M C8 )sin(C 9 + C lf ,M) (A6) 
'max c 

Fits for two specific airfoil sections are shown by the dashed lines in 
figure A. For the V23010 - 1.58 (type 1 above), only the first three terms 
of equation (A6) were required. However, for the VR-7 (type 2 above), all 
ten coefficients are nonzero. In both cases, the functions were obtained by 
use of a least-squares curve-fit program. 

Reynolds number effects- Figure 5, from reference 2, presents the com- 
bined effects of M and for two airfoil sections. It can be seen that the 
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of fort of increasing Rjj Is to shift t ho complete 


i’j vs M I'lirvr now, Mils, 
'max 

with tho curves tending to collapse together above about Rfj * I'lO 1 ’ In ono 
case, an<l ahovo about Rfj ** b'llV' for tho othor. Therefore, tho oonstant l’j 
In oqu.it Ion (Ah) ahovo will ho determined (from low M oxpor'montal .lata for 
tho particular soot Ion 1 as a function of tho oxpootoil R^ range. 

Negative angle of aitaok- Kxpor imont al ilata lor o. at negative angle 

" max 

of attaok, as a funot Ion ol M, aro not avallahlo for most alrfo'l soot Ions. 

It Is assumed that tho variation of c .> with M at negative a has tho 

max 

same form as for positive Tho sign of o. Is negative, of ootirso, and 

mux 

for oamborod sections tho oonstant term C ( In equation (Ah) will have a dif- 
ferent magnitude as wol 1 as a different sign. The it j term Is obtained trom 
tho experimental lift curve at low M. If no data aro available tor negative 
i\ oven for low M, It can bo estimated as 


max 

<*, ) 


[(•■« ) 

- * (■'*’ . ) , 
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That is, for a cumbered section, negative c ■ Is lower In magnitude than 

max 

the positive c ( hv an amount double the lift Increment due to the can, her 

max 

The accuracy ol this estimate was evaluated hv comparing it with measured 
values for several airfoil sections, as presented In table S. The error Is 
less that 10T In six of the seven cases studied. 

l.AKUK ANf.I.K OK ATTACK 

Very little experiment.*' section lift coefllcient data for angle ol 
attack more than a few degrees above the stall angle are available. Sect Ion 
lilt coefficients through '180° for the NAt’A 0012 and the NAOA MAO 12 ate 
presented in figure h. A curve fit ol these data in four segments Is given 


bv 

- 1.1 - 1 . 78|0. 01745(a) - 0. 785 1 1 ’ , 2 1' . a - d()“ (A8) 

c v - -l.l + 1.78[0.01745(a) - 2. lSf> ) •' , 0(1“ . a - lbO“ (A**) 

c v , - -0.7M , 160“ v „ v 172.5“ (A10) 

Cj - - 5 . 82 [ n - 0.01745(a)) , 172.5“ . a L 180“ (AH) 


Kquat Ions (A8 ) through (All) are taken from reference 7, except that (A8) Is 
begun at a - 22“ rather than lb” since this yields better agreement with tint 
experimental data in figure b. 

It is assumed that airfoil section lift coefficient at large angles of 
attack is an odd. symmetric function about a » 0 (even for cambered sections). 
Thus e v for large negative angles is given bv the above equations, except 
that the sign of c is revers'd: 
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,’ V-l 


i-_ -aa.xa?*- 


1—— 1-— T-X 1 J — L.J 


i ..J_ -X ... 


■U — l L i- J 



- -1.1 + 1 . 78 1 0. 01 7AS ] 

1 - 0. 78S *1 , 

-1) 1° V - J J ° 

( A 1 2 ) 


- i.i - i. ; hi o. oi 7 vs t , ! 

. - 2. ISb) ' . 

-lb('' v a .. -40" 

(A 1 ») 

'' V' 

- 0. 7b 1 , -172.S" - 

a _ -lb0° 


(A 1 t ) 


- S . 8 2 | t - 0. 1)1 7 s Vj a | ) 

. -180“ _ 

- 1 7 2 . r > ° 

( A 1 r > ) 


F.qual Ions (AH) through (A1S) art 1 .is .mim'd to hold lor .ill airfoil soot Ions, 
regardless of camber, thickness, M, or Kjg. 

lhie to tlio nn-ancr amount of tost data for o ,• .i'hht stall, tho soot Ion 
lltl oooltlolont is assumod to bo a straight lino botwoon o .• and 

i • 22“. Tlion, from oijuat Ion (AH), m ' lx 


0.81 » 


0.81 3 +■ 


•> > _ 


(22 - a) 


stall 


stall 


(Alb) 


whorl* It Is assumod that c .. occurs at 

'max 


V 

max 

‘stall " do 7 /da + 


(Al?) 


If tho particular airfoil section exhibits gradual stall ch.uactcrlst los, 

equation ( A 1 7 ) does not accurately predict tho angle at which c occurs. 

max 

However, equation (A17) may still bo used In equation (Alb), because tho pur- 
pose of tho latter Is to approximate the lift curve in the region between tho 
(assumed) linear lift-curve region and the high angle of attack region 
(a j : 22 °). For negative angles of attack. 


'' “ -°- 8n + 22 ♦ a 


(c, ) + 0.81 

\ max/- , 


- (22 + a) . 




nog stall 


a v (A18) 

nog stall 


where 

('■ ) 

\ max/-., 

’nog stall dc j/da + “d (A19) 
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Al’VKNDIX H 


1>K.R 1VAT ION OK UKAi: 1'OK.KK U‘ l K.NT EOPATIONS 


SMAU. ANOl.K OK ATTAt'K 


Kor airfoil sect loo angle ot <Ut.uk less t turn the stall .mu i *' . section 
ilru*!. eoeff lelent Is essent Li 1 1 v constant wttli Mach number below t lie drag 
divergence Mac li number Mp|). Above Mpp, * tie l *‘ 'R itses verv steeple. llms 
It Is assumed that 

'.! * ^Vm-h * M - M m> * '* ~ 'stall (H ' 

t . 

'j • "aV, * a« ;M ' V ■ ” V- ' '.till 

where (e .) , vie. AIM. and M are t Hint Ions ot sect Ion angle >'t attack. 

I* m~ o \i nn 

!>!;!>; i\u K t t I *' ifitl at tow M 

Section drag coct t le t i at tor smooth alrtolls at low Midi mimhei . with 
tullv turbulent houndai v Liver (the important ease tor most practical .ipi* l l 
eat ions), iiuiv he estimated iii't , l) hv 


is 1 (: , ' V • 

A ‘ \ v) "i <,;»'»* 


Fi| nation t it * ) was derived lor symmetrical alrtoll sections; howovet , as shown 
hv experimental data piesentcd in reterence camher has a small cited on 
the minimum drag coot l le lent . Further, an Inspection o t the section data in 
reterence I stiows that minimum oj occurs at approx Imnt o l v > * 0 *ot hath 
svmiiH't r lea 1 and eamhered sect ions. l'heietote, <is discussed In reteience 1, 
the i j vs i relat L'u obtained tor a symmetrical alrtoll mav he applied with 
good accuracy to an alrtoll with the same thickness distribution but a cam- 
bered mean line. The factors S.^, l./e, e^/cp, aitil k are eonslants ti'i a 
sped tie alrtoll section, and are obtained t‘ rom graphs provided In leler- 
ouce 1. Those graphs are reprodm'ed here <is tigures 1 through U*. I’he turbu- 
lent skin t fiction ilt .ig coefficient Is represented hv 

0.4SS 


('" K K \.,,) ' 


reterence 10. The effective K [yj (ref. 1) is given hv 


K Mp(5 !•) ,Sa • 


i i i 
• «» 





win* tv Kfj , is based on fm'-striM# voloeltx and j i i r l . i 1 eliotd, oxrrpl wnero 
it Is spec i I Ifni 1 v taken as BMP 1 ' (‘bird term ot oq . (Hi)!. 


Pr.ig Oivorp.ouee Ma. li Numbei 

Tlx dray. tl i vomonre Marh number of an airfoil soot loti is dot I nod as t lit* 
Marti number lor whirl) (dr ( |/dM) " 0.1 as airspeed i . increased at ran . taut 
sort I on anyle ol allark. Measured dray J i vr rcenee Marh nuirbeis for sovrtal 
airfoil sort Ions arc presented In t lyuro II. lbe Mpp data points in 
flyure 11 haw horn fitted with straight linrs, or romh i na t ions of two straight 
linos tn t lit* c.imbnod airfoil eases. Corre'.Mi <n Is r«*rv pood lot t hr NAi'A 
4-dlylt, r *-diyit, hl-series, and Wort maun airfoils, t l.ronphout tin* t ranr.r 
for whirl) data woo av.iilahlo. However, t ho data to*' tin* N At ' A htA-serles 
airfoils have a shift to different st ro 1 yht - 1 lm seymeni s ahovr ibout anyle 

of attack. Also, t In* NAI'A b4A-sor los data appoar to hav»* an ahrnpt tluelua- 
t Ion at about a “ 10". It was not ost ah 1 i shod whrthrr a yeoer.il frond of 
6-sor'os airfoils Is represented hv thoso data ( 1 1 so noto that flu* NAi’A h‘> > 


data are well represent* 

*d bv a st r 

a lyht 

1 illo) . si lire* 

a yeneral express 1 on was 

desired, it was assumed 

Hint M W) 

ran 

bn represented 

bv oqn. 1 ions of tin* form 


M PP “ A 

+ Hi 

. *» - 

PI) 

( Bh ) 


M nn “ ( 

4- |)a 

. ‘ - "|, 

PP 



where iip,... is tlu* pork of tho Mpp data. For svmmrt r Ira 1 sort Ions, 
is zero, C » A, and I) “ -B. Studv of oxpor Imont a I data in relereneos 1 and • 
indiratos that Mpp Is novor loss than 0. 1, reyard 1 ess ot anyle ot aft. irk. 


Slopo of r j Curve Above* M ^ 

Tho slopo of tho i‘j rurvo abovo Mpp was measured from oxpor imont a ! 
data for several airfoil soitlons, and tho rosults aro prosotitod in I Ip.nro 1-’. 
The curves plot tod in tho fiyure wore obtaiiiod from 1 oast -squares rurvo 1 its 
of tin* oxpor Imont a l data. Th.* 4-diylt sorlos airiolls aro easilv roprosontod 
by i ttrvos ot tho form 

Jc d 

■ A + Ba + CiV 
un 

llowovor. tho trend of (doj/dMV Is more complex for tho other two soot Ions 
studied. Curve fits of tho form 

dr 

~ - A + Ba + Ctr' + Pa * 
dM 

wort* roqulrrd to dupl Irate tho rise and then level injt-o ft of the experimental 
data. Note that the quadratic equations will turn upwards, and the ruble 


16 



LARGE ANGLE OF ATTACK 


Experimental see lion drag coefficients through *180° are presented for 
two NACA .lirtoll seetlons tn figure 1 l. These data are companion to the lift 
coefficient data In figure h. The curve fit of the data shown In figure \\ 

Is given bv • 1 

c d - M8(;sln a|)'- . IS* i | a I _ 180° (BUI 

The lower limit of IS* for equation ( B1 2 > is arbitrary, but reflects the fact 
that the cl tecta ot section tvpe and of Mach number are significant onlv for 
snull a. Equation (BUI was obtained from reference l, hut the constant 
factor has been increased slightly to provide better correlation with the 
experimental data in figure 1\ Equation (BUI is assumed to he applicable to 
all airfoil sections, regardless of section thickness, Reynolds number or 
Mach number. 

It' section angle of attack is greater than a stj|11 , hut less than IS®, 
the drag coefficient is assumed to he a straight line between and 

0U19 - (e , 

- 1IA l * 

St 


d 'lV - 


'statl 

i (la - 

all ! 


stall 
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APPENDIX C 
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DERIVATION OF MOMENT COEFFICIENT EQUATIONS 
SMALL ANCLE OF ATTACK 


The variation of airfoil section pitching moment (about the quarter-chord) 
with Mach number Is small below the moment divergence Mach number, and mav be 
neglected In that region. Also the moment is essentially Independent of 
Reynolds number. For most airfoils, the slope of the moment coefficient curve 
is approximately constant (often zero) with angle of attack until 1 the stall 
angle is reached, whereupon c m breaks sharplv. Therefore, it i9 assumed 
that 


dc 


d 


m , 

~ (u) , 


M s*s® 


< n > - "stall 


(Cl) 


where c m is c m at a • 0. At positive stall, c breaks in the negative 
(nose-down) c m direction. At negative stall, the reverse is true. 


Moment Divergence Macn Number 

An Inspection of experimental airfoil section pitching moment coefficient 
plotted against Mach number (e.g., ref. 4) shows that, at a fixed angle of 
attack, c m is essentially constant with M until a certain M Is reached, 
whereupon c m diverges rapidly with further increases in M. A pitching 
moment divergence Mach number is defined such that M Mn is the Mach number 
for which ! dc m /dM =0.5 as airspeed is increased at constant section angle 
of attack. Using this definition, M^ as a function of angle of attack was 
obtained from experimental data for several airfoil sections, and is plotted 
in figure 14. Although there is some scatter, the data are reasonably well 
fitted by straight lines, or combinations of straight lines. Therefore, 

Mj^p will be represented by equations of the form 

Mj^p - A + Ba , a > a p (C2) 


Mm * C + Da , a < a 

MD 


(C3) 


where u Pj^p Is the peak of the data. 

C « A, and D * -B Experimental data 
Mjfl) is always greater than about 0.3, 


For symmetrical sections a P - 0, 

*MD 

in references 3 and 4 indicate that 
for any angle of attack. 
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Above Stall or 

There are verv few data available for pitching moment above stall or 
Mjqp, and no general trends are discernible. It is assumed that ’ m varies 
linearly with a between n s| . ( j| anti 20° (assumed to be the start ot the 
large u region). A linear relation of e TO with M is also assumed between 
!%) and M = 0.95. At M = 0.95 it Is assumed that o m has roaehe.1 t he 
same large- i- region value of -0.077 as at < = 20°. The selection of >1 = 0.9t 
for this relation Is arbitrary. Thas, 

(cm) + 0.077 

c - -0.077 + (70 - a) , 20° " u " * , (04) 

m 20 - a . . st .i l I 

Still 1 




< v "'hd + °- 077 


0.95 




(M ‘ m md ) 


m md 


( 05 ) 


Note that the condition a n s t a 11 , M Miqp may occur. In which case equa- 
tion (04) is used to evaluate (e m )|,|p in equation (05). For negative angles 
of attack, these equations become 


0.077 


(c m>. 


0.077 


st ]1 1 


20 + a 


(20 + a) , 


20 ° 


i 


stall 


Vail 


(06) 


and 


(c m*MD 


< c m>„D - °- 077 


°- 95 ’Mmd (M " ^ 


M ' % 


(07) 


where a , , is negative, 
stall 

LARC.E ANCLE OF ATTACK 


Experimental section pitching moment coefficients through *180° are pre- 
sented in figure 15 for the NACA 61A012 and NACA 0012 airfoils. A curve 
through the data was determined in four sections, as shown in figure 15. The 


curve is given by 

c m - -0.00802 (a - 20) - 0.077 , 20° < a < 67° (08) 
c m - -0.619{sln(C. 0260a - 1.26)1 0 - 3 ^ , 67° < a <. 162° (C9) 
C;,, - -0. 00838 (a - 162) - 0.320 , 162 6 < a < 170° (CIO) 
c m - 0.0387 (a - 170) - 0.387 , 170° < a < 180° (CJ1) 


The lower limit of 20° for these large-angle equations is arbitrary. 

20 
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4. 
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Section pitching moment at large angle of attack is assumed to be an odd, 
symmetric function about a * 0, even for cambered sections. Thus c tor 
large negative ang io •* (s given bv the above equation*, except that the sign of 


c m is reversed (c m is positive) 

c m - O.uti802(j .( 20) t- o.o?; , - a 7 0 a -20° (i 12 ) 

c m - 0. 61 9 1 s- in (0. 0260 1 a | - 1.76) J°- v » 8 , -162° a - -67" (CM) 

c m - 0. 008 >8 ( !o j - 162) + 0. 320 . -170° < a < -162° (CU) 

c m - -0. 0387 ( | a j - 170) + 0.387 , -180° v a < -170' (CM) 
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AiTKNnix n 

OOMITTFR PROCRAM FOR FVAl.lATION OK A! RIOT I. 
SFCT ION AKROPYNAM IC CliAKACTFK 1 ST ICS 


A digital computer program vas written to a 1 1 ow rapid cv.ilu.it T>u at' t ha 
airfoil sort loti fora*’ and moment roof I' ie i ant aquations prasantad In this 
report. The program calculates seat Ion litt, drag, and pitching moment ! - 
flcients for angles of attaak from -ISO 1 ' to + 180 ! ' l and for a range of Much 
numbers (M 1.0) which rn.iv ha selected bv the user. Input pai imat rrs for the 
specific airfoil typo under studv are required. ..sample input lor the 
NACA 23012 airfoil st>ation is presented in table 3. t*ala ulalod output is in 
tlie form of lift coof fie lent , drag aoeffiaient, and pitching moment coeffi- 
cient tables. In each table, the calculated eoefl icieiit is (H inted is a luna- 
tion of section angle o' attack and of Mach number. Calcul ited aerodynamic 
coefficients fi'r the NACA 2 101.' are plotted In figures 3 anl I. These data 
were calculated by the p’s gram as a rosin of the Input data shewn in t »M«* 3 . 
A FORTRAN listing of the computer program is presented in tnis appendix. 
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-1 -l-J ..3- 


Input data 


TABLE I . - DATA REQUIR EMEN TS 


(dc ,, /dl '*) inc 

t/c 


M,, 


( l * 0 ) i . 


c e 

coef f ic ients 

max 


C 1 

» C ](1 for +o case 

r ! 

■* Cj 0 for -a case 


I E quat ion n umbe rs 

Lift coefficient 


(A2), (A3), also f 1 gs 
(A3), (A4 ) , (A5) 

(A4), (AS) 


Drag coefficient 


(B6), (B7 ) 
(B6) , (B7 ) 
(B5) 


M dd coefficients A, B. C, D 

Op 
* nn 


dc d /dM coefficients A, B, C, D (B8) 

°c (B8) , (B9>, (Bll) 

(BIO), (Bll) 


Moment coefficient 


coefficients: A, B, C, D |(C2), (C3) 

(C2) , (C3 






TABLE 2.- SAMPLE INPUT FOR NACA 21012 


Symbol 

FORTRAN 

name 

Value 

1 — 

Symbol 

FORTRAN 

name 

Value 

— 

NMNP 

10 

L/c 

RLC 

2.035 

— 

SMACH 

0- 

S A 

SSA 

1.18 

— 

DMACI1 

.1 

c s /c F 

CSCF 

.037 

t/'c 

TC 

.12 

K 

FK 

1.55 

R N 0 

RNO 

8.10*10 f> 

Coefficients for M DD : 

(dc j/da) lnc 

CLA1 

. 100/deg 

A 

AMDD 

0.730 

M 1 

AMACH1 

.80 

B 

BMDD 

-.0246 

(«(,>! 

ALFZl 

-1.20 deg 

c 

CMDD 

.830 

M r 

AMACH2 

IT 

GO 

D 

DMDD 

.0246 

(a 0 ) 2 

ALFZ2 

-. 70 deg 

a p 

ALFPDD 

-2.0° 


r DD 



Coefficients for +c <? 

J 

Cl 

ALFCUT 

10.0° 


'max 

c 

Cl 

Cl 

1.622 

a R 

ALFREF 

-2.0° 

c 2 

C2 

. 337 

Coefficients for 

dc d /dM: 

c 3 

C3 

-2.316 

A 

AC DM 

0.274 

c„ 

C4 

.0 

B 

BCDM 

.0253 

C 5 

C5 

.0 

C 

CCDM 

.00273 

C & 

C6 

.0 

D 

DC DM 

.000264 

Ci 

C7 

.0 

Coefficients for 

m md : 

C« 

C8 

.0 

A 

AMMD 

0.810 

Cg 

C9 

.0 

B 

BMMD 

-.026 

c 10 

CIO 

.0 

C 

CMMD 

.910 

Coef f ieient s 

for -ci > : 

C ma*c 

D 

DMMD 

.026 

Cl 

C1N 

-1.200 

ap MD 

ALFPMD 

-2.0° 

C 2 

C2N 

-.250 

c *o 

CMO 

-.010 

c 3 

C3N 

1.716 

dc m /da 

DCMDA 

.0014/deg 

c<« 

CAN 

.0 




C 5 

C5N 

.0 




C 6 

C6N 

.0 




Cl 

C7N 

.0 




c B 

C8N 

.0 




c 9 

C9N 

.0 




ClO 

CION 

.0 
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TABLF 1. 


MEASURED ANP 
FOR SEVERAL 


b 


Reference 2. 
Equation (A2). 


CALCULATED 
AIRFOIL 


LIFT-CURVF .M.Oi’F 


Airfoil 

Much 

Measured* 1 

Calculatt 

sect Ion 

number 

dc t /d i 

dc,/d.» 




per dof^ 

por iio> 

NACA 0012 

0. 30 

0.103 


i 

1 

.40 

. 108 

0.108 

1 

j 

.60 

.128 

.127 

! 

.80 

! . 200 

.197, 

j NACA 23015 

. 30 

. lOu 



.40 

. 105 

.105 


.60 

.120 

.125 


.75 

.155 

.170 

j NACA 6 j AO 12 

. 30 

.100 

___ 

! 

. 39 

. 106 

.104 

i 

.58 

.125 

.121 

) 

1 

.74 

.160 

.157 

| VFRTOI, 

.30 

. 122 


V23010-1. 58 

.40 

. 124 

. 129 


.60 

.145 

.151 

1 

. 77 

.240 

.206 

VF.RTOL VR-7 

. 30 

.110 

___ 


.40 

.118 

.115 


.62 

.138 

.139 



.75 

■ ,80 __l 

.178 


1)1 f f erenc 

•y 


0.0 

.8 

-3.0 

.0 

4.2 

9.7 

-1.9 

-3.2 

-1.9 

4.0 

4.1 
•14.2 

-2.5 

.7 

- 1.1 
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TABLE 4.- MEASURED MACH NUMBER EFFECT ON ATRFOIL 
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71 1 I r 1 7~f~ H 1 1 . n~i :i.-.L-L.L,m 
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20 *0 


60 80 

«, deg 


!40 160 180 


Figure 


1.2 ULl 
-20 0 


100 12° 


i.- 


C ale vi 


1 »ted lift. dra>\, a,lJ 
NM'A 23012 airtoll 


p l teh i nH _n,t ' raeMt 
section. «t M 


coeff to lent* 

0 . 1 . 


tor t'u' 


3 b 


1 


.1. u J H XIX J 




M = 0.4 r 


ODA MEASURED, REF. 2 
CALCULATED 


0 8 



A 


.04 



OD A MEASURED. REF 2 
— CALCULATED 









OOA&Q MEASURED, REF. 4 



(c) Section pitching. moment coefficient. 
Figure 3.- Concluded. 
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MACH NUMBER 



(a) Airfoil section V23010-1.58. 
(b) Airfoil section VR-7. 
Figure 4.- Maximum lift coefficient. 
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MEAN VA LUE O F AIRFOIL 
PRESSURE COEFFICIENT. S 



PERIMETE 



i Figure 8.- Airfoil section perimeter/chord ratio. 

i 

J 



5^ .08 

o 

< .04 


Figure 9. 





M 00 M FOR WHICH dCy/dM • 0.1, AT CONSTANT n 




N'Ai'A f.vnw»>tr U'.il fiuir-J 1*1 1 series .ilrloll sections. 

<b> NAi'A e.imbered iowr-Jiglt series U t rf of l sections. 

Ki^urf II.- HimsuioiI strfol 1 sect Ion Jr.«»; divergence Much number. 


mjjn 



13 V J_] 1.1 IJJuJJ 






• 

1 

.8 


M PO 


.6 


D0 = M FOR WHICH dC d /dM - 0.1. AT CONSTANT u 


O NACA 64A010 . REF 12 
A NACA 64A310 , REF. 12 


O NACA 65 2 -215, £ = 0.5. REF. 11 
AWORTMANN FX69-h-098 , REF. 4 



(e) NACA 64A series airfoil sections. 

(f) NACA 65 series and Wortmann airfoil sections. 
Figure 11.- Concluded. 
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rv— U. '“'i -J L .. J 1... J i . .1 I J 




dC d /dW = 0.274 + 0.0253 « + 0.0C273 a 2 - 0.00264 a 3 

(7 X 



0.466 - 0.0261 a + 0.00137 <v 2 
0.399 -0.0378 « + 0.00245 a 2 


0.147 + 0.0127 a + 0.01 1 1 a 2 - 0.000874 Q 3 


O NACA 0012 34 , REF. 3 
A NACA 231 2, REF. 3 
□ NACA 23012, REF. 4 
[7 FX69-H-098, REF. 4 


a; deg 


8 


10 


12 


Slope of 


l ’d curve Above drug divergence Mach number. 
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